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Abstract. The nature of the pairing state in iron-based superconductors is the 
subject of much debate. Here we argue that in one material, the stoichiometric 
iron pnictide KFe 2 As 2 , there is overwhelming evidence for a d-wave pairing state, 
characterized by symmetry-imposed vertical line nodes in the superconducting gap. 
This evidence is reviewed, with a focus on thermal conductivity and the strong impact 
of impurity scattering on the critical temperature T c . We then compare KFe2As 2 to 
Bao.6Ko.4Fe2As2, obtained by Ba substitution, where the pairing symmetry is s-wave 
and the T c is ten times higher. The transition from d-wscve to s-wave within the same 
crystal structure provides a rare opportunity to investigate the connection between 
band structure and pairing mechanism. We also compare KFe 2 As2 to the nodal iron- 
based superconductor LaFePO, for which the pairing symmetry is probably not rf-wave, 
but more likely s-wave with accidental line nodes. 
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(a) Ba .6Ko.4Fe 2 As2 (b) KFe 2 As 2 




Figure 1. Sketch of the Fermi surface of the iron pnictide Bai_ a; K :c Fe2As2 in the 
k z = plane, for two K concentrations, (a) For x = 0.4: two hole-like cylinders at 
the zone center (r) and two electron-like cylinders at the zone corner (X)[l]. (b) For 
x = 1.0: three hole-like cylinders at the zone center (r), one hole- like cylinder near 
the zone corner (X), and no electron pocket [2, 3]. Recent ARPES data [4] confirms 
this Fermi surface topology in the k z — plane; it also reveals the presence of a small 
3D pocket centered at the Z point of the Brillouin zone (at k z = ir/c, not shown), 
consistent with calculations [2]. The disappearance of the electron pocket in going 
from x = 0.4 to x = 1.0 is thought to play a major role in the dramatic change in the 
superconducting state: from s-wave symmetry with T c ~ 40 K to d-wave symmetry 
with T c ~ 4 K. 

1. Introduction 

Iron pnictides show that high-temperature superconductivity can be achieved in a good 
metal on the border of antiferromagnetic order. A key question is the role of the 
antiferromagnetic spin fluctuations in causing the pairing. Early on it was proposed 
that the two-band nature of the materials (Fig. 1), with one (or more) hole-like Fermi 
surface at the zone center (r point) and an electron-like Fermi surface at the zone corner 
(X point), is responsible for both the spin-stripe order and the superconductivity [5]. 
Pairing would rely predominantly on the inter-band interaction that links the hole and 
electron Fermi surfaces, producing a pairing state with s-wave symmetry but with a gap 
that has opposite sign on the two Fermi surfaces, the so-called s± state. 

In this Article we consider a particular iron pnictide, KFe2As2, which offers an 
unprecedented window on the workings of magnetically-mediated superconductivity, as 
it undergoes a phase transition from one pairing symmetry to another, induced simply 
by substitution of Ba on the K site. The change from ci-wave symmetry in KFe2As2 to 
s-wave symmetry in Bai_ a; K E Fe2As2 is accompanied by the appearance of the electron 
pocket in the Fermi surface and a ten- fold increase in T c (see Fig. 2). This change 
of symmetry was predicted theoretically by renormalization-group calculations [6] and 
recently confirmed by experiment [7]. Here we review the evidence. 
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Figure 2. Phase diagram of the iron-pnictide superconductor Bai_ a; K :c Fe2As2 
as a function of K concentration x. At x = 0, the parent compound is 
antiferromagnetic (AF) below a temperature Tn which decreases with x (black 
squares [8]). Superconductivity exists in a dome-like region below the transition 
temperature T c . With increasing x, T c rises to reach a maximal value of 38 K at 
optimal doping [x ~ 0.4), outside the region of coexisting antiferromagnetic order 
(shaded in grey). At x = 1.0, single crystals have the lowest disorder of any iron-based 
superconductor. The properties of KFe2As2 are very well characterized, and T c = 4.0 
K in the cleanest samples. In this article, we compare the system at two different 
concentrations: x = 0.4 (blue dot) and x = 1.0 (red dot). We show that the symmetry 
of the superconducting state changes from s-wave at x = 0.4 to d-wave at x = 1.0. 
This means that there must be a phase transition in between. 



1.1. d-wave vs s-wave 

There are two kinds of nodes where a superconducting gap goes to zero: those that 
are imposed by the symmetry of the pairing state, and those that are accidental, not 
imposed by any symmetry but the result of an anisotropic pairing interaction. In a 
d-wave state, of B lg symmetry in a tetragonal structure, the gap goes to zero along 
diagonals in the x-y plane. In the simplest 2D model, with A = Aocos20, A(0) = 
when the azimuthal angle <fi equals 7r/4, 37r/4, 5ir/4 and 7tt/A. The order parameter 
changes sign as goes from to 7r/2 and so it breaks the tetragonal symmetry of the 
lattice. If the Fermi surface is a quasi-2D cylinder located at the center of the Brillouin 
zone, then the gap will necessarily have 4 vertical line nodes running parallel to the axis 
of the cylinder. See Fig. 3. 

An s-wave gap has A\ g symmetry and so preserves the rotational symmetry of the 
underlying lattice. It can have accidental nodes, but the number of nodes where the gap 
changes sign must preserve tetragonal symmetry. On a quasi-2D Fermi surface cylinder 
in a tetragonal structure, an s-wave gap with nodes (called an "extended s-wave" gap) 
can go to zero along lines that can either be vertical or horizontal. For example, the 
gap can have 4 pairs of vertical line nodes, with 4 negative minima separating 4 positive 
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Figure 3. Top panels: Sketch of the superconducting gap A((p) (red dashed line) 
as a function of azimuthal angle in the basal plane of a tetragonal quasi-2D 
superconductor with a circular Fermi surface (blue solid line). A d-wave gap (on 
the left) changes sign 4 times as it goes around the Fermi surface and as a result the 
gap necessarily goes to zero at four nodes (small black circles). It breaks the four- fold 
rotational symmetry of the lattice. An s-wave gap (right) does not break the four-fold 
rotational symmetry of the lattice. In general it does not change sign, but it can do so, 
at accidental nodes, as sketched here for example. Bottom panels: Effect of increasing 
the impurity scattering rate T. For a d-wave gap (left), the nodes always remain, as 
they are imposed by the symmetry of the order parameter. Impurity scattering rapidly 
suppresses the gap magnitude. For an s-wave gap (right), scattering makes the gap 
more isotropic and will in general remove the nodes. It is less effective in reducing the 
gap magnitude. 

maxima, as sketched in Fig. 3. 

2. Effect of impurity scattering on T c 

Impurity scattering has a qualitatively different effect on <i-wave and extended s-wave 
states, as pointed out early on in the study of cuprate superconductors [9]. Scattering 
cannot remove the <i-wave nodes, as they are imposed by symmetry. By contrast, 
scattering in the s-wave state will in general make the gap more isotropic, and eventually 
remove (lift) the nodes (Fig. 3). 

Intra-band non-magnetic scattering is very effective in suppressing the T c of a d- 
wave superconductor. The scattering rate V needed to suppress T c to zero is of the order 
of k-BT C Q/h, where T c0 is the clean-limit value of T c . In a single-band model with unitary 
scattering, the critical scattering rate r c is given by HT C = 0.88 k^T^ [10, 11]. Intra-band 
scattering is much less effective in suppressing the T c of an s-wave superconductor [12]. 

In this context, let us examine the behaviour of iron-based superconductors. In 
Fig. 4, we plot T c vs T, both normalized by T c0 , for various 122 superconductors. In 
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Figure 4. Effect of impurity scattering on the transition temperature T c of various 
iron-pnictide superconductors. The scattering rate T is estimated using the residual 
resistivity pQ. Data for BaFe 2 As2 (squares), SrFe2As2 (diamonds) and CaFe2As2 
(triangles) optimally doped with Co, Pt or Ru are reproduced from ref. [17]. Data 
for KFe2As2 (red circles) are reproduced from ref. [7]. Both axes are normalized by 
T c0 , the value of T c in the clean limit (r 0); T c0 = 3.95 K for KFe 2 As 2 [7] and 
T c0 = 26 K for optimally doped BaFe 2 As 2 , SrFe 2 As 2 and CaFe 2 As 2 [17]. All lines 
are the same theoretical Abrikosov-Gorkov curve for the decrease of T c with scattering 
rate, plotted as T c /T c0 vs r/r c , where the critical scattering rate T c for suppressing 
T c to zero is the only adjustable parameter. For KFe2As2, hT c = 1.3 /cb?cO ( see 
text), while for the other compounds, TiT c = 45 fee^co- The 30- fold difference in 
T c /T c0 is compelling evidence for superconducting states of different symmetry. Inset 
: Zoom on the KFe2As2 data at small T. The black dashed line is the theoretical 
expectation for a d-wave superconductor, assuming a single-band Fermi surface, for 
which h,r c = 0.88 k B T cQ [10, 11]. 

KFe2As2, KT C ~ k^Tco, in agreement with expectation for a ci-wave state. T c is estimated 
using the critical value of the residual resistivity po f° r which T c will go to zero, namely 
^cnt ^ 4 Q cm j n XFe2As2 [7]. Using the normal-state specific heat coefficient 
7n = 85 ± 10 mJ/K 2 mol [14, 15] and the average Fermi velocity v p — 4 x 10 4 m/s [13], 
we obtain TiY c = h/2r c ~ 1.3 ± 0.2 ksT C Q, via the normal-state thermal conductivity 
« N /r = L /p c rit = 7 N^r c /3, where L = (tt 2 /3)(e/k B ) 2 . 

Note that in this and other estimates below, we use single-band models and relations 
to compare theory and experiment in KFe2As2, a multi-band material. The reason why 
a single-band approximation may work reasonably well, is that the Fermi surface sheets 
that dominate the spectral weight are the three relatively similar hole-like cylinders 
centered on the F point (Fig. 1) [2, 13]. The weak magneto-resistance of KFe2As2 was 
offered as evidence of single-band behavior [16]. 

In dramatic contrast, T c is suppressed 30 times more slowly in the materials based 
on BaFe 2 As 2 , SrFe 2 As 2 and CaFe 2 As 2 , optimally doped by substitution on the Fe site 
(with Co, Co/Cu, Pt, Ru, Ni, and Pd). Indeed, all of those were found to fall on the same 
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curve, with T c going to zero at TiT c ~ 45 /cb^co [17]. This is completely incompatible 
with rf-wave pairing, and indeed some of these materials have gaps without nodes, e.g. 
Co-doped BaFe2As2 [18, 19]. This 30-fold difference between KFe2As2 and materials 
that otherwise have an identical crystal structure and a very similar Fermi surface is 
compelling evidence that the pairing symmetry must be different. We now use heat 
transport at very low temperature to show that the pairing symmetry of KFe 2 As2 is 
indeed o?-wave. 

3. Theory of heat transport in a d-wave superconductor 

The thermal conductivity k measured at very low temperature is a powerful probe of 
the superconducting state [20]. It is a bulk measurement, insensitive to modifications 
that might occur at the surface of a sample. It is directional, and so can access the 
anisotropy of the gap structure. And it probably has the best energy resolution of all 
experimental probes for distinguishing nodes in the gap from deep minima, dictated by 
the lowest temperature of the measurement, typically 50 mK or so. Finally, as we shall 
see, it can distinguish between symmetry-imposed line nodes and accidental line nodes. 

The presence of a residual linear term in k(T), whereby k/T is non-zero in the 
limit of T — > 0, is unambiguous evidence of zero-energy fermionic quasiparticles. In a 
metal or in the normal state of a superconductor induced for example by applying a 
magnetic field greater than H C 2, this residual linear term, labelled kh/T, is given by the 
Wiedemann-Franz law : 

£ = £ ■ <u 

T po 

In the superconducting state at zero magnetic field (H — 0), a non-zero residual 
linear term, labelled kq/T, is unambiguous evidence that the superconducting gap must 
vanish somewhere on the Fermi surface. The gap must have nodes, whether points or 
lines, accidental or symmetry-imposed. By contrast, if the gap does not have nodes, 
even if it might have deep minima, then k,q/T = 0, as confirmed by measurements in 
the s-wave superconductors V 3 Si [21], MgB 2 [22], NbSe 2 [23] and LuNi 2 B 2 C [24], for 
example. 

The theory of heat transport in a rf-wave superconductor is well developed [10, 
25, 26, 27, 28, 29]. The zero-energy quasiparticles associated with the nodes in the 
d-wave gap produce fermionic heat transport in the T — > limit. These quasiparticles 
are induced by impurity scattering, and further excited by a magnetic field or a rise in 
temperature, as we shall now describe. 

3.1. Magnitude of k /T and its dependence on impurity scattering 

Calculations show that the residual linear term in a <i-wave superconductor is 
independent of scattering rate and phase shift [25] , and free of Fermi-liquid and vertex 
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corrections [29]. For a quasi-2D superconductor in the clean limit (hT <C A ), it is given 
by [25, 29]: 

^00 _ *| 1 (2) 

T 3h c va 

where c is the interlayer separation, vp is the in-plane Fermi velocity, and v a is the 
quasiparticle velocity tangential to the Fermi surface at the node, proportional to the 
slope of the gap at the node. We will use this expression to estimate the theoretically 
expected value of k /T for the cuprate superconductor YBa 2 Cu 3 07. 

Eq. 2 may also be expressed in terms of the residual linear term in the normal-state 
electronic specific heat, 7n, and the slope of the gap at the node, /xA (expressed in 
terms of the gap maximum, A ) [25]: 

Koo h 7n^f f o\ 
T ~ vr ^A • [6) 

In the simple case where A = A o cos(20), /j, — 2. Assuming the weak-coupling value of 
the gap, A = 2.14 k B T c o, and \x = 2, we get: 

^oo = h 7n^| 

T ~ 4.28vr k B T c0 ' [ ' 

We will use this expression to estimate the theoretically expected value of kq/T for 
KFe 2 As 2 . 

The remarkable feature of ci-wave superconductors (or indeed any superconductor 
with symmetry- imposed line nodes) is that Kq/T is independent of the impurity 
scattering rate V. This comes from a compensation between the impurity-induced 
growth in quasiparticle density and the proportional decrease in mean free path. This is 
what is meant by "universal heat conduction" , first verified experimentally in the cuprate 
superconductor YBa 2 Cu 3 07 [30]. Although it is strictly true only in the clean limit 
(%T <C A ), the change in kq/T is expected to remain weak up to hT/k B T c0 ~ 0.5 [10]. 

3.2. Dependence of k/T on temperature 

So far, we have discussed the limit T — )■ and H — )■ 0, where nodal quasiparticles 
are excited only by the pair-breaking effect of impurities. Raising the temperature 
will further excite nodal quasiparticles. Calculations for a <i-wave superconductor show 
that the leading-order finite-temperature correction to Eq. 2 (or Eq. 3) is a positive 
contribution to the electronic thermal conductivity that grows as T 2 [25, 26]: 

T ~ {1 + a ^f 

where a is a dimensionless number and 7 is the impurity bandwidth, which grows with 
the scattering rate T [25]. In the limit of unitary scattering, 7 2 oc T; away from unitary 
scattering, towards Born scattering, •y 2 grows more rapidly with T [25]. 



+ , (5) 
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3.3. Dependence of kq/T on magnetic field 

Increasing the magnetic field is another way to excite quasiparticles in the 
superconducting state. In the absence of nodes, the quasiparticle states at T = 
are localized in the vortex core and heat transport occurs only as a result of tunneling 
between adjacent vortices [20]. The dependence of Kq/T on H is then exponentially 
slow. If the gap has nodes, however, the field will shift the energy of delocalized states 
outside the vortex cores and enhance the residual density of states, which will grow 
typically as y/H. This causes an immediate rise in Kq/T with applied field [20, 27, 28]. 
Calculations for a quasi-2D d-wave superconductor in the clean limit (HF <C A ) reveal a 
non-monotonic increase of Kq/T vs H, with most of the increase to reach k^/T occurring 
in the last quarter of H/H c2 [28] (see Fig. 9). 



4. Experimental data for KFe 2 As 2 

In this section, we compare the properties of KFe 2 As2 to those expected of a d- 
wave superconductor, both from theory, as outlined in the previous section, and from 
observation in cuprates, the archetypal ci-wave superconductors. We also compare to 
theoretical calculations for an extended s-wave gap applied to pnictides [12]. Finally, 
we compare with the equivalent data on Bao.6K .4Fe2As2. 

We begin by showing, in Fig. 5, thermal conductivity data for KFe 2 As 2 and 
Ba .6K .4Fe 2 As 2 , at H ~ 0. (A very small field H = 0.05 T is applied to kill the 
superconductivity of the contacts [19].) It is clear by inspection that there is a large 
residual linear term in KFe 2 As 2 and none in Bao.6Ko.4Fe 2 As 2 . The fact that kq/T = 
in the latter [31] is sufficient proof to rule out ci-wave symmetry for the superconducting 
state at x = 0.4. 



4-1. Magnitude of Ko/T 

When extrapolated to zero field (see inset of Fig. 9), the residual linear term in KFe 2 As 2 
is [7]: 

Y = 3.6 ± 0.5 mW/K 2 cm . (6) 

Let us compare this with the value expected for a d-wave gap. Using Eq. 4 and 
7nU 2 = 2.3 ± 0.3 x 10 13 mJ/K 2 cm s 2 [7], we get 

?Y = 3.3 ± 0.5 mW/K 2 cm , (7) 

in excellent agreement with experiment. 

Let us also compare with experimental data on cuprate superconductors. 
Measurements on YBa 2 Cu 3 07 near optimal doping, at a doping p ~ 0.18 where 
T c = 89 K give [32]: 

^ = 0.16 ± 0.02 mW/K 2 cm . (8) 
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Figure 5. Temperature dependence of the in-plane thermal conductivity n of 
Bai_ a; K a ;Fe2As2, measured in a very small magnetic field (H ~ 0), plotted as k/T 
vs T 2 , for two K concentrations: x — 0.4 (blue) [31] and x — 1.0 (red) [7]. The lines 
are a linear fit, used to extrapolate the residual linear term at T = 0, labelled kq/T. 

Using Eq. 2 and the ratio of quasiparticle velocities measured by ARPES in 
Ba2Sr 2 CaCu208+<5, vf/va = 16 ± 1 at p ~ 0.16 [33], yields 

1° = 0.16 ± 0.01 mW/K 2 cm , (9) 

since the interlayer separation c = 5.85 A in YB^CusOy. This is in perfect agreement 
with experiment. The same is true for Tl2Ba2Cu06+<5 at p = 0.19 with T c = 84 K: the 
measured value, kq/T = 0.08 ± 0.02 mWK 2 cm [34], is half that in YBCO because the 
interlayer separation, c = 11.6 A, is twice what it is in YBCO. 

It is important to appreciate, however, that Eqs. 2 and 3 are only valid in the 
clean limit, when hT <C A . This condition holds well for the cleanest samples of 
KFe2As2, with RRR values of 1000 or so [7]. An easy way to quickly estimate the 
ratio hT/A is via the ratio of superconducting-state to normal-state residual linear 
terms: (k /T)/ (kn/T) ~ hT/A . For a sample of KF 2 As 2 with p = 0.2 [ifl cm, 
this gives TiY / Aq = 0.025. The clean-limit condition almost certainly also holds for the 
cuprates YBa 2 Cu 3 0j / and Tl 2 Ba 2 Cu0 6+( 5 near optimal doping, but it doesn't in strongly 
overdoped Tl 2 Ba2Cu06+5 [35]. Neither does it hold for most other cuprates, including 
Ba2Sr2CaCu20 8 +«5 [36] and Lai_ x Sr x Cu0 4 [21, 36] (see discussion in [37]). In the dirty 
limit, the thermal conductivity is no longer a simple measure of the o?-wave gap near 
the node, but instead it reflects the normal-state character of the conductivity. 

These quantitative comparisons with theory show how reliable the absolute values 
of Ko/T can be in assessing whether a particular superconductor has a o?-wave gap. There 
are two reasons for this. First, Ho/T does not depend on the concentration or type of 
impurities or defects in the sample. Secondly, there is typically a direct connection 
between the slope of the gap at the node, v&, which controls the density of excited 
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Figure 6. Dependence of kq/T and on impurity scattering rate T, normalized by T c q, 
the disorder-free superconducting temperature. Data for the pnictide KFe 2 As2 (red 
circles; from [7]) and the cuprate YBa 2 Cu3C>7 (blue squares; from [30]) are normalized 
by the theoretically expected value for a c?-wave superconductor, kqq/T = 3.3 and 0.16 
mW/K 2 cm, respectively (see text). The typical dependence expected of an s-wave 
state with accidental nodes is also shown, from a calculation applied to pnictides (black 
line; from [12]). 

zero-energy quasiparticles, and the gap maximum A , and hence T c . Neither of these 
reasons applies to the case of accidental line nodes in an extended s-wave gap. Indeed, 
while Eq. 2 still holds [12], v& is now strongly dependent on the impurity scattering 
and its value is not related in any simple way to T c . So the remarkable agreement 
in KFe2As2 between the measured no/T and the value estimated via Eq. 4 would be 
entirely fortuitous in a nodal s-wave state. 

4-2. Dependence of kq/T on impurity scattering 

A powerful test of whether line nodes are imposed by symmetry, as in a <i-wave 
superconductor, or accidental, as in an extended s-wave state, is the dependence of 
kq/T on impurity (or defect) scattering. In Fig. 6, we see that kq/T in KFe2As2 is 
unchanged by a 10-fold increase in T (i.e. a 10-fold increase in p ) [7]. The classic data 
on the d-wave superconductor YBa2Cu30y [30], reproduced here for comparison, show 
a similar insensitivity to disorder. More than any other measurement, this rules out 
the possibility of an accidental line node in KFe2As2. Calculations based on a two-band 
model for the pnictides produce a wide variety of curves for Kq/T vs T [12]; a typical one 
is reproduced in Fig. 6. The value in the clean limit is different from the ci-wave k Q0 /T, 
by an amount which depends on parameters. In general, kq/T will depend strongly on 
T. Here k,q/T starts by growing, as t>A falls with increasing T, then it goes through a 
near-divergence as v a — > just before the nodes lift and finally it plummets to zero 
when the nodes are eventually lifted by sufficiently strong scattering (see Fig. 6). 
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It is probably not possible, and certainly artificial, to adjust the various parameters 
of a two-band (or N-band) model of s-wave superconductivity in the pnictides to ensure 
that the following four facts are simultaneously reproduced: 1) correct value of Kq/T; 
2) constancy of kq/T vs scattering up to T/T c = 0.5; 3) suppression of T c to zero at T c ; 
and 4) correct value of r c (i.e. ~ T c0 ). As we saw, these four properties are naturally 
obtained in a <i-wave state, with no adjustable parameter. 

4-3. Dependence of k/T on temperature 

Another test of d-wave theory applied to KFe 2 As2 is the temperature dependence of 
k/T at low temperature (T <C T c ). As seen in Fig. 5, the data for KFe 2 As 2 obey 

^ = f(l + AT 2 ) , (10) 

with A = 23 K' 2 . 

Before we attribute this AT 2 term to electronic excitations, we need to know what 
is the phonon contribution. Indeed, while phonons do not contribute to the residual 
linear term, they do of course carry heat, and k p oc T a , with a typically between 2 and 
3 [21, 37]. There are several ways to show that the phonon term k p is at least 10 times 
smaller than the AT 2 term below 0.4 K, and so cannot account for it. 

The maximal value of the phonon conductivity k p at sub-Kelvin temperatures 
is achieved when the phonon mean free path reaches the size of the sample. If 
that boundary scattering is diffuse (on rough surfaces), the T dependence will be 
k p = BT 3 , or Kp/T = BT 2 [37]. Sutherland et al. have calculated the expected phonon 
term k p from the known phonon properties of the iron-based superconductor LaFePO 
(specific heat and velocity) and the given sample dimensions [46]. Since their sample 
dimensions are comparable to ours and the sound velocities of iron-based materials are 
comparable, their estimate also applies approximately to our KFe 2 As2 sample. They 
obtain B — 1.0 mW/K 4 cm, while our measured slope is 100 mW/K 4 cm. There is 
no way that such a huge slope could be due to phonons. It must therefore be due to 
electronic excitations. 

A second way to rule out phonons as the cause of the T 2 slope is to compare 
with a sample that has 10 times more impurity scattering. This should not affect k p 
at those temperatures, where the phonon mean free path is controlled by the sample 
boundaries [37]. In Fig. 7, we see that it does. Indeed, the slope of k/T drops by at least 
a factor 10. A third way is to look at the data on Bao.6K .4Fe2As2, in Fig. 5. Because 
Ba .6K .4Fe2As2 has no nodes, electrons do not contribute to the slope of k/T either 
(for T < T c /100). We see that k p ~ BT 2 , and B is 10 times smaller than the slope A 
in KFe 2 As2- 

Having shown that the AT 2 term in KFe 2 As 2 is electronic, and hence due to nodal 
quasiparticles, we have yet another, independent confirmation of <i-wave behavior. As we 
saw, according to Eq. 5, a T 2 dependence is the expected low-temperature behavior. A 
T 2 dependence was indeed observed in YBa 2 Cu 3 7 , with k/T = {k /T){1 + 19 T 2 ) [32]. 
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Figure 7. Thermal conductivity of KFe2As2 for a clean sample (closed symbols), 
whose po increased by a factor 1.15 over a period of 6 months (before, red circles [7]; 
after, blue squares), and a dirty sample (open black circles [41]), whose p is 10 times 
higher than that of the clean sample. The data is plotted as k/T vs T 2 , normalized by 
the value (ko/T) in the T = limit. Note that the actual value of kq/T is identical in 
all three data sets, within error bars (at H — > 0). These data show that the T 2 slope 
is strongly reduced by an increase in impurity or defect scattering, as expected of a 
d-wave superconductor (see Eq. 5). 

This is yet another difference with respect to an extended s-wave state, for which 
calculations do not obtain a quadratic growth in k/T at low temperature [12]. 

Moreover, in <i-wave theory the coefficient A is expected to decrease with scattering, 
as we just saw it does. In the limit of unitary scattering, A oc 1/T , so that a 10-times 
larger V would yield a 10-times smaller slope [25], as roughly observed. We can be 
more precise by using data on a high-quality sample whose RRR dropped slightly over 
a 6-month period, by a factor of 1.15. The data is shown in Fig. 7. The increased 
disorder scattering did not change k /T within error bars, but it did change the T 2 
slope noticeably, by a factor 1.60. So we see again, now within the same sample, that 
an increase in V causes a drop in the T 2 slope while keeping k /T constant, as expected 
of a d-wave superconductor. 

The temperature below which the T 2 dependence of k/T in KFe2As2 sets in, 
T ~ 0.1 T c , is a measure of the impurity bandwidth 7. It is in excellent agreement with 
the temperature below which the penetration depth X a (T) of KFe2As2 (in a sample with 
similar RRR) deviates from its linear T dependence [2]. This is beautifully consistent 
with the theory of c?-wave superconductors [38]. 

4-4- Dependence of kq/T on magnetic field 

The presence of nodes in the gap will cause Kq/T to be immediately enhanced by 
application of a magnetic field H. A rapid initial rise in kq/T vs H has been observed 
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Figure 8. Residual linear term, kq/T, normalized by the normal-state thermal 
conductivity in the T = limit, Kn/T, as a function of magnetic field H applied 
along the c axis, normalized by the upper critical held H C 2, for two concentrations of 
Bai_ x K x Fe2As2: x — 0.4 (closed blue circles) [31] and x = 1.0 (open red circles) [41]. 
The values of H C 2 used here are 100 T and 0.8 T, respectively. Data on KFe2As2 are 
from a sample in the dirty limit [41]. For comparison, we reproduce corresponding 
data for the isotropic s-wave superconductor Nb (open grey squares [20]) and the dirty 
d-wave superconductor Tl-2201 (open grey circles [35]). 

in the cuprate d-wave superconductors YBa 2 Cu 3 07 [39] and Tl 2 Ba 2 Cu0 6+( 5 [35], for 
example. In the dirty limit, measurements on Tl 2 Ba2Cu06+5 (with T c = 13 K) [35], 
reveal a smooth monotonic increase of kq/T vs H all the way to H c2 ~ 7 T, as reproduced 
in Fig. 8. In dirty-limit samples, KFe 2 As2 [41] and Tl2Ba2Cu06+<5 [35] show very similar 
curves of k.o/T vs H/H c2 (Fig. 8). Note that the similar value of (kq/T)/(k^/T) ~ 1/3 
at H = reveals a similar level of scattering, with HT/A ~ 1/3. 

In the clean limit, one expects no/T to fall abruptly as the field is reduced below 
H c2 . This is true for any type of gap. It comes from the fact that vortices enter 
the sample and immediately introduce an extra scattering process, limiting the long 
quasiparticle mean free path of the normal state. In Fig. 8, this behavior is seen in 
the classic s-wave superconductor Nb. A behavior almost identical to that of Nb was 
observed in clean-limit samples of the iron pnictide LiFeAs [40] , showing this particular 
stoichiometric pnictide to have an s-wave gap. (At H = 0, it also has Ho/T = 0.) 

Fig. 9 shows how theoretical calculations for a clean d-wave superconductor contain 
these two features: immediate initial rise at low H and sharp drop just below H c2 [28]. 
Data on a clean sample of KFe2As2 behave precisely in this way, and are in remarkably 
good agreement with <i-wave calculations (see Fig. 9). In the case of cuprates, a full 
comparison with <i-wave theory has not so far been possible in the clean limit. Although 
the clean limit is certainly reached in YBa2Cu3C>7, measurements have been limited to 
H < H c2 . 

In Fig. 8, the data for Bao.6K .4Fe2As2 are also shown, and the increase in k /T at 
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Figure 9. Field dependence of Ko/T measured in two samples of KFe2As2 for a 
current parallel to the a axis (J||a) and a field parallel to the c axis (if ||c): one in the 
clean limit (r/r c = 0.05; full red circles [7]), the other in the dirty limit (r/r c = 0.5; 
open red circles [41]). The black dashed line is a theoretical calculation for a d-wave 
superconductor in the clean limit (Sr/Ao = 0.1) [28]. Inset: Zoom of the same data 
at low field (left axis). Also shown is kq/T measured for a current along the c axis 
(J || c) in a sample with T/r c = 0.10 (black squares [7], right axis). 

low field is seen to be very flat [31]. This is further evidence that there are no nodes 
in the gap, and it also shows that there are no deep minima. Because the same flat 
dependence is seen for a heat current along the c axis, it reveals a full isotropic gap at 
x = 0.4 [31]. 

4-5. Dependence of k /T on current direction 

Our assignment of a <i-wave gap in KFe2As2 must withstand yet another stringent test: 
the anisotropy of the gap structure. As sketched in Fig. 3, a <i-wave gap on a single 
quasi-2D cylindrical Fermi surface (at the zone center) would necessarily have 4 line 
nodes that run "vertically" along the length of the cylinder (along the c axis). In such 
a nodal structure, zero-energy nodal quasiparticles will conduct heat not only in the 
plane, but also along the c axis, by an amount proportional to the c-axis dispersion of 
the Fermi surface. In the simplest case, c-axis conduction will be smaller than a-axis 
conduction by a factor equal to the mass tensor anisotropy (t>§ in Eq. 3). In other 
words, (kcq/T)/ '(K a o/T) ~ (K c -^/T)/(K a ^/T) = (o- c /o- a )^, the anisotropy in the normal- 
state thermal and electrical conductivities. 

This is confirmed by calculations for a quasi-2D d-wave superconductor [42] , whose 
vertical line nodes yield an anisotropy ratio in the superconducting state very similar 
to that of the normal state. 

This is what we see in KFe2As2 (inset of Fig. 9): (n a o/T) / '(kco/T) = 20 ± 4, very 
close to the intrinsic normal-state anisotropy {o~ a /a c )^ = 25 ± 1 [7]. This shows that the 
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Figure 10. Sketch of different scenarios for line nodes on the Fermi surface of 
KFe2As2 (x — 1). In panels (a) and (b), the effect of a line node on the anisotropy 
of heat transport is illustrated in a single-band model, whereby the Fermi surface 
is approximated by a single warped cylinder. The c-axis dispersion (along k z ) is 
responsible for conduction along the c axis, producing an a-c anisotropy in the normal- 
state (N) conductivity a, (<7 C / <7 a )N, which in KFe2As2 is equal to 0.04 [7]. In the normal 
state at T — >■ 0, the very same anisotropy is observed in the thermal conductivity k, 
so that (k c /kci)n = 0.04. The onset of superconductivity will alter this anisotropy, 
depending on the anisotropy of the gap. For a gap with a (horizontal) line node 
in the basal plane (a), thermal conduction in the basal plane will be much better 
than conduction along the c axis, so that in the superconducting state (S) in zero 
field, (K c /K a )s <C (k c /k q )n- By contrast, if the line node is vertical, lying within 
the a-c plane, then the anisotropy of the normal state is more or less preserved, and 
( K c/ K a)s — ( k c/ k <i)n [42]. In (c), we illustrate how in a two-band scenario, a vertical 
line node present only on one of the two Fermi surfaces can also lead to a much reduced 
anisotropy in the superconducting state. This will happen if c-axis conduction in the 
normal state is caused predominantly by that surface (with strong 3D character) which 
does not have the accidental line node. 

nodes in KFe2As2 are vertical lines running along the c axis. 

If there were only a single quasi-2D sheet in the Fermi surface of KFe2As2, the 
anisotropy could equally be due to symmetry-imposed vertical lines in a ri-wave gap or 
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accidental vertical line nodes in an extended s-wave gap. However, the fact that the 
Fermi surface of KFe2As2 contains several sheets with very different c-axis dispersions 
provides compelling evidence in favor of rf-symmetry. The argument is as follows. 
As sketched in Fig. 10c, if one Fermi surface sheet dominates the c-axis transport 
in the normal state because it has the strongest c-axis dispersion, but it does not 
have nodes in the superconducting state, then the anisotropy (« / k c )s measured in 
the superconducting state will be much larger than the anisotropy (« / k c )n in the 
normal state. Now the Fermi surface of KFe 2 As 2 does contain a small pocket at the 
Z point whose strong 3D character [2, 4] will make it a dominant contributor to c- 
axis conduction. In an extended s-wave scenario, there would be no particular reason 
why the gap would develop nodes on this small pocket, and so the anisotropy would 
be very different in the superconducting and normal states, unlike what is measured. 
By contrast, in ci-wave symmetry this pocket would automatically and necessarily have 
nodes, as any other zone-centered sheet, thereby ensuring that transport anisotropy 
remains similar in the superconducting and normal states. The same argument applies 
to the three T-centered cylinders, one of which has a much stronger c-axis dispersion 
than the other two [2, 4]. Again, all three must have vertical line nodes to preserve the 
similarity of anisotropies in superconducting and normal states. This is automatically 
ensured by <i-wave symmetry, but is not the case in s-wave scenarios such as that of 
ref. [43], where only one of the three sheets has nodes. 

5. LaFePO: a nodal s-wave superconductor 

In light of what we know of KFe 2 As 2 , it is interesting to investigate the case of the 
iron-based superconductor LaFePO. In some ways, LaFePO shares several similarities 
with KFe 2 As 2 : it is stoichiometric, with a relatively small T c ~ 7.5 K, and its gap has 
nodes, as detected by a linear T dependence of the penetration depth [44] and a residual 
linear term in the zero- field thermal conductivity [45] . Are those line nodes imposed by 
symmetry or accidental? Is the symmetry d-wave or s-wave? 

These questions were addressed in a recent study of heat transport in LaFePO [46], 
where some of the aspects considered here were covered. Two features are singled out 
by the authors as inconsistent with rf-wave symmetry. First, T c seems to be robust 
against impurity scattering. Although only one sample was studied, the estimated 
scattering rate for that sample is such that hT ~ A; B T C . <i-wave superconductivity would 
be almost entirely suppressed by such a large scattering rate. It doesn't appear to be. 
This robustness certainly contrasts with the high sensitivity of KFe 2 As 2 , and is more 
reminiscent of optimally-doped BaFe 2 As 2 , SrFe 2 As 2 and CaFe 2 As 2 [17]. Note that a 
systematic study of the impact of impurity scattering on the T c of LaFePO is needed to 
confirm this claim. 

The second feature is the temperature dependence of the electronic part of k/T, 
which is linear at low temperature (T < T c /15) [46] rather than T 2 as in KFe 2 As 2 and 
YBa 2 Cu 3 7 . These two facts suggest that an extended s-wave scenario is more likely 
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for LaFePO, in agreement with renormalization-group calculations [47, 48]. A good test 
would be to study the dependence of k,o/T on T and on current direction. 

Note that the iron pnictides BaFe2(Asi_ :r P :r )2 [49] and Ba(Fei_ x Ru x )2As2 [50] also 
have nodes in the gap. The robustness of T c to impurity scattering in the latter [17] is 
again an argument in favor of a nodal s-wave gap. 

6. Summary 

In this section, we summarize the evidence presented in this Article that leads us 
to identify the pairing symmetry of the stoichiometric iron pnictide superconductor 
KFe2As2 as d-wave, and why the pairing symmetry of Ba-doped KFe 2 As2 must be 
different. It is impressive that the respective d-wave and s-wave symmetries of these 
two otherwise quite similar materials was correctly predicted by renormalization-group 
calculations [6]. 

6.1. BaosKoAFe2As2: an s-wave superconductor 

The fact that k /T = at H = in Bao.6Ko.4Fe 2 As 2 for current directions parallel 
and perpendicular to the c axis [31] shows that there are no nodes in the gap of this 
material, immediately ruling out d-wave pairing. The very slow increase of k /T with 
field, again in both current directions [31], also reveals that the gap has no deep minima 
and is essentially isotropic. This suggests s-wave symmetry, with or without a change 
of sign between the hole and electron Fermi surfaces. The striking fact about this 
superconductor remains its high critical temperature, T c ~ 40 K, and its very high 
upper critical field, H c2 ~ 100 T [51, 52]. 

6.2. KFe 2 As2-' a d-wave superconductor 

The six independent experimental properties that confirm KFe2As2 as a d-wave 
superconductor are listed in Table I. We summarize them here. T c is suppressed by 
impurity scattering at a rate consistent with %Y c /k-BT C Q ~ 1.0. This rapid rate, precisely 
as expected for d-wave, is 30 to 50 times more rapid than what is observed in other 
iron-pnictide superconductors with the same crystal structure. This is a huge difference, 
pointing compellingly to a difference in pairing symmetry. Nodal quasiparticles have 
a conductivity which is unaffected by a 10-fold increase in impurity scattering. This 
universal heat conduction is the signature of line nodes imposed by symmetry. The 
magnitude of that conductivity and its a-c anisotropy are precisely those expected of 
a d-wave gap, with vertical line nodes on all the zone-centered Fermi surfaces. The 
initial growth in k/T with temperature is precisely T 2 , precisely as expected of a d- 
wave superconductor. Moreover, the T 2 term is rapidly suppressed by scattering, again 
as expected. The entire field dependence of k /T, from H = to H = H c2 , is in striking 
agreement with calculations for a d-wave gap. 
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Table 1. Quantitative and qualitative characteristics of superconductivity in 
KFe2As2, compared to the corresponding characteristics expected of a <i-wave 
superconductor (with symmetry-imposed vertical line nodes) and of an extended s- 
wave superconductor (with accidental vertical line nodes). T c is the critical impurity 
scattering rate needed to suppress T c to zero. T c0 is the clean-limit value of the 
superconducting transition temperature. kq/T is the residual linear term in the 
thermal conductivity extrapolated to T = and H = 0. kqq/T = h-f^v^/(4:.28TrkBT c o). 
{n c /n a ) s is the anisotropy of heat conduction in the superconducting state at T — > 0, 
namely the ratio of kq/T measured along the c axis over kq/T measured along the a 
axis. (kc/k<i) n is the corresponding anisotropy of heat conduction in the normal state, 
equal to (cr c /<7 ) N , the anisotropy in the electrical conductivity a. 
* Extended s-wave with accidental vertical line nodes, 
f Value obtained in optimally-doped ^4i?Fe2As2 (see Fig. 4 and [17]). 
| Strongly dependent on details - can be much more or much less than 1.0. 





KFe 2 As 2 


d-wave 


s-wave* 


hTj k B T c0 


1.3 ±0.2 


0.88 


45 f 


(k /T) / (k o/T) 


1.1 ±0.5 


1.0 


t 


{K c / K a ) s / (k c / K a ) N 


1.25 


~ 1.0 


t 


T dependence of k /T 


weak 


weak 


strong 


T dependence of k/T 


T 2 


T 2 


T 


H dependence of Kq/T 


rapid 


rapid 


rapid 



By contrast, an s-wave gap with accidental line nodes on one Fermi surface cannot 
account for all six of these features simultaneously. Impurity scattering will in general 
have a strong effect on Kq/T and it is known to have a weak effect on T c in s-wave 
superconductors of the 122 structure - the opposite of what is seen in KFe2As2- The 
magnitude of Hq/T in clean samples will only be right by accident, and only by another 
coincidence will the a-c anisotropy be right. The T dependence will not in general be 
the observed T 2 . 

We may now ask how the change in order parameter symmetry from s-wave at 
x = 0.4 to d-wave at x — 1.0 proceeds as a function of x in Bai_ x K x Fe2As2. One 
scenario is two separate superconducting phases, separated by a non-superconducting 
region. A second scenario would be to have a continuous superconducting phase where 
the change in symmetry from s-wave to ci-wave proceeds via an intermediate state of 
(s ± id) symmetry, a complex order parameter that breaks time-reversal symmetry [53] . 
Further work is required to elucidate this question. 
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Note added in proof - By adding Co impurities in KFe 2 As 2 , a recent study [54] 
has confirmed that T c falls rapidly to zero with impurity scattering, roughly at p cr u = 
4.5 fj.fl cm, and the residual linear term in the thermal conductivity is indeed universal, 
remaining approximately constant even when the normal-state conductivity is decreased 
by a factor 30. 
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